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Abstract 

We analyze a stimulated revival (echo) effect for the breathing modes of the 
atomic oscillations in optical lattices. The effect arises from the dephasing 
due to the weak anharmonicity being partly reversed in time by means of 
additional parametric excitation of the optical lattice. The shape of the echo 
response is obtained by numerically simulating the equation of motion for the 
atoms with subsequent averaging over the thermal initial conditions. A qual- 
itative analysis of the phenomenon shows that the suggested echo mechanism 
combines the features of both spin and phonon echoes. 
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I. INTRODUCTION. 



The coherent manipulation of the atomic center of mass motion in optical lattices |T]] 
by means of the nonstationary off-resonant dipole potentials has provided new experimen- 
tal capabilities to study dynamical systems with time-dependent potentials [0- 0. In the 
off-resonant coherent regime, the parametric non-adiabatic excitations of the optical lattice 
give rise to oscillations of the atomic momentum and coordinate distribution dispersions 
(breathing modes) and may be used for the manipulation of the coordinate or momen- 
tum dispersions of the atomic distribution by means of squeezing in phase space [Q. 

The squeezed states of matter in various systems have been extensively discussed in 
connection with the possibility of overcoming the standard quantum limit for noise imposed 
by vacuum fluctuations (see e.g. 0,0] )• In an optical lattice, the squeezing effect can be 
significant in the classical regime producing classical squeezed states p| and resulting in the 
reduction of thermal fluctuations. 

Recently, the breathing modes of atomic oscillations in the optical lattices have been 
observed experimentally with the use of the Bragg scattering techniques based on the fact 
that the cross section of the Bragg diffraction scales with the Debye- Waller factor which is 
related to the dispersion of the atomic coordinate distribution . The observed decay of 
the oscillations may be due to both the dissipative and the dephasing effects. 

In classical oscillating systems with anharmonicity the frequencies have a continuous 
distribution due to the energy-dependent corrections. This leads to the dephasing of small 
oscillations P] . Together with possible irreversible dissipation, the dephasing determines an 
apparent damping of the oscillations. The dephasing, caused by the nonlinearity effects, 
can lead to partial revivals under certain circumstances. Examples of such revivals are the 
phonon (see in 0), spin and photon echo effects [|10||, |[TT[| . 

One has to distinguish between the classical stimulated revival (echo) caused by external 
perturbations applied to the oscillating system, and the quantum revival related to the 
discrete structure of the energy levels of the system which does not require any external 
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disturbances [§,0. As will be shown below, these two revivals can have different time 
scales and therefore can be separated. 

In this paper, we show that an analog of the phonon echo effect exists for the breathing 
modes of atomic oscillations in the optical lattices. However, the proposed echo mechanism 
possesses new features compared with the conventional phonon echo mechanism. First, the 
proposed mechanism is related to the oscillations of the dispersion of the atomic distribution 
as opposed to the average coordinate and momentum oscillations involved in the phonon echo 
model [^. In contrast with the conventional model of phonon echo [^, the excitation of the 
system in our model is achieved by means of the parametric modulation of the frequency 
of the atomic oscillations as opposed to applying a time-dependent external force to the 
atoms 1^. As shown below, the suggested echo mechanism incorporates the features of both 
spin (photon) |lTO| , p!I| and the phonon echo mechanisms due to the specific nature of the 



parametric excitation of the oscillations. This effect can be employed to study the revival 
and damping mechanisms of the breathing modes in optical lattices. 

II. THE MODEL AND NUMERICAL RESULTS. 

An atom subjected to a off- resonant laser field experiences an energy shift of the ground 
state proportional to the intensity of the field. If the field is formed by a standing wave with 
large detuning, the effective potential for the atoms in the ground state is given by |ll3 



U{x) =V[l-cos{2qx)], (1) 

where x denotes the atomic center of mass coordinate, V is the amplitude of the dipole 
potential (proportional to the intensity of the laser field), and q is the wave vector of the 
laser field. Normally, the initial state of the system is formed by Doppler cooling, so that the 
trapped atoms are in thermal equilibrium near the minima of the optical potential and form 
an optical lattice. In this paper we will consider the case of temperatures T much larger 
than the energy of the atomic oscillations at the bottom of the potential. Consequently, no 
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quantum effects will be taken into account. On the other hand, the temperature should be 
low enough to insure that the atoms do not escape from the potential well. If the magnitude 
of the optical lattice potential is time-dependent, then 

V = Vo[l + k{t)], (2) 

where k{t) describes the parametric excitation of the lattice (Vo=const). Therefore, the 
effective potential can be approximated as an oscillator with weak anharmonicity and time- 
dependent frequency. The harmonic frequency corresponding to the potential (|l|) is given 

by 

Uo = 2qJ—, (3) 
V m 

where m is the atomic mass. We assume that the parametric excitation of the system occurs 
in the form of two short pulses 

k{t) = SiS{t-ti) + S2S{t-t2), (4) 

where Si_2 and ti^2 characterize the intensities and the instants when the pulses are applied 
to the optical lattice. In what follows, we assume that ti = and t2 = tp. The equation of 
motion is given by 

£ = -i sin(2£)(l + 265(r) + 265(r - r^)), (5) 

where we introduced the dimensionless variables 

x(r) = gx(r), r = ujot, (6) 

with ^1 = uJoSi/2, ^2 = ^082/2, Ti = ujQti = and T2 = Tp = Uotp. 

The physics of the echo effect in the optical lattice can be understood as follows. We 
assume that initially the atomic system is in thermal equilibrium characterized by the values 
of coordinate and momentum dispersions mujQ < >t=< >t /m = ksT. The first 
pulse initiates the squeezing of the atomic phase space distribution creating the classical 



squeezed states 0. In harmonic systems, such squeezing results in oscillations of < a;(t)^ > 
and < p{t)'^ > with the frequency 2ujq. The non-linearity leads to dephasing and therefore 
to the apparent damping of these oscillations. Since the dephasing is reversible in time, 
one expects to observe at least a partial revival of the oscillations after the application of 
the second pulse to the system. A qualitative physical picture of the phenomenon will be 
discussed in the next section. Here we concentrate on the numerical results. 

In order to study the echo effect, we performed a numerical simulation using the equation 
of motion (|^) with subsequent averaging of x^(r) over the initial conditions x(0), x{0) 
distributed in accordance with the initial thermal ensemble. The results for the normalized 
average dispersion < x^(r) > / < >t are presented in Fig.l. One can see the expected 
damping of the squeezing oscillations after each pulse due to the dephasing effect. The 
characteristic damping rate is increasing with an increase of the temperature. The most 
remarkable feature of the time evolution of < x^(r) > is the pronounced echo effect at 
T ^ 2Tp. The revival time is approximately 2rp according to the general considerations 
involved in the two-pulse phonon and the spin echo models. Note however, that the revival 
of < x'^{t) > is only partial (approximately 0.25 of the initial amplitude). 

In Fig. 2, we present the dependence of the echo amplitude on the intensity of the second 
pulse. The echo amplitude is very small for small values of ^2, then it increases sharply in 
the intermediate region of .^2 and saturates in the nonlinear regime. 

III. QUALITATIVE ANALYSIS. COMPARISON WITH THE PHONON AND 

SPIN ECHO. 

In order to clarify the mechanism of the proposed echo effect and compare with the other 
echo mechanisms, we consider the approximate analytical solution of Eq.(^ valid for the 
low temperature (small anharmonicity) . In this case, it is sufficient to take into account 
only the first (quartic) correction to the effective harmonic potential. Equation (^) reduces 
to 
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X = -x{l + 2ei<5(r) + 265(r - t,)) + - i^ (7) 

where we neglected the modulation of the non-linear term assuming that the amplitudes of 
the pulses and ^2 are sufficiently small. For r 7^ 0, Tp, the Eq.(^ has an integral of motion 
corresponding to the energy of the atomic system E = E/AVq (i.e., there are 3 different 
values of energy E^, = 0, 1, 2 for times r < 0, < r < Tp and r > r^, respectively). 

For small anharmonicity and amplitudes of the external pulses, the solution of the Eq. (|^ 
for 7(r) = x^(r) up to the second order terms in E]^ can be obtained in the form 

i{t) = Ck + [Afcexp(2iu;(^fc)r) + 5^ exp(4it^(Efe)r) + c.c], 

(8) 

Ck = Ek + jEl, 

where k = 0,1,2. From Eq.(7), it follows that = Ek{l + exp(2z(/)fc)/2 and 
Bk = —Al/6, where (pk corresponds to the phase of the atomic oscillations in different time 
intervals. Note that Ak are the complex constants characterizing both the amplitude and 
phase of the atomic oscillations and uj{E) stands for the effective frequency. This frequency 
is determined in the limit i^^ — as 

ujiE) = 1-Ie. (9) 

The solution given by Eqs.® and (P) corresponds to the first term of the asymptotic ex- 
pansion for the anharmonic oscillator in the limit of weak anharmonicity and large evolution 
time [|14| . 



Each pulse in Eq.(^ results in a jump of 7 and E, while 7 remains continuous. We 
obtain from Eq.(|^ 

7(r+) - 7(t-7) = -4a7(r.); E{r^) - E{t;) = 24%(r.) - ^7(^7), (10) 

where = 0, Tp and and r~ denote the moment of time just before and just after each 
pulse and k = 1,2. After the application of the first pulse, we obtain from Eqs.(||) and (|TUp 
for small ^1 
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Ai=Ao + iUEa + 2Ao) (11) 

and 

Ei = Eq-2iUAq-AI). (12) 

Combining Eqs. (|TTD , (|^), and averaging over the random initial phase and the thermally 
distributed initial energy in the limit of low temperatures (0 = T/4Vo << 1), we obtain the 
time evolution of the coordinate dispersion in the form 

<^(r)>=e + e[z67^^^^^^ + c.c.], rp>r>0, (13) 

with the dephasing rate of oscillations defined by 

= e. (14) 



Eq.(13) indicates that the amplitude of the oscillations of the coordinate dispersion exhibits 
an apparent damping. As discussed above, this damping is caused by the anharmonicity 
which results in the dependence of the effective oscillation frequency on energy, and the 
thermal energy distribution leading to an averaging out of the oscillations. Therefore, this 
damping is reversible in time. The dephasing rate given by Eq. (|14D is in a good qualitative 
agreement with the numerical results discussed above. Fitting the data presented in Fig.l 
for < r < Tp by Eq. (|T^) gives the value t^^ ~ 1.10 for the dephasing rate in our numerical 
simulations. 

Now we will analyze the effect of the second pulse. Applying Eq.(p!OD for = and 
Ts = Tp consecutively, one obtains the amplitude A2 in Eq.(|D after the second pulse (t > Tp)m 
the limit of small ^2 

A2{Tp) ^ ii(rp)(l + 2^6) + ^6^1 - li^2EiAl{Tp) (15) 

where 74i^2(t) = ^1^2©^*'^^^^'^''^ correspond to the envelope solution for 7(r) given by Eq.(^). 
The y4- vectors rotate with the frequencies 2u!{Ei^2) in the complex plane. The energy after 
the second pulse is given by 



E2 = E,- 2ze2(Aie2^"(^i)-- - c.c), (16) 

with Ai and Ei defined by Eqs.(|ll]) and (|12D. 

In order to obtain a qualitative picture of tlie eclio effect, we disregard the energy change 
due to both pulses upon the frequencies lj{Ei^2)- Eq. (p!5|) indicates that the amplitude A2{Tp) 
after the second pulse is a linear combination of the amplitude before the pulse ~ ^i(Tp) and 
the amplitude A2e{Tp) = —i^2EiAl{Tp)/3. Note, that A2e{Tp) is proportional to the complex 
conjugate of Ai{Tp). 

In Fig. 3, we present the two vector amplitudes Ai and A'l rotating with the frequencies 
2uji > 2uj'i in a complex plane. If the two vectors have the same phase = at r = then 
the phase difference between Ai and A'^ is accumulated over the time interval Tp. The phases 
of the complex conjugate vectors A\ and A( are inverted with respect to Ai and A'^. The 
vector A'-^ was rotating faster than the vector Ai and therefore A'* is delayed with respect 
to A\. If the vectors continue rotating with the same frequencies 2uJi and 2uj[, the phase 
difference between the vectors Ai and A\ increases contributing to the dephasing effect. In 
contrast, the phase difference between A\ and A( vanishes at the time r = 2Tp giving rise 



to the echo analogous to the well known spin and photon echo effect |10|, [11|. Taking into 

account Eqs.(|8|)- ([T0|) , ([T5|) one arrives at the echo contribution after thermal averaging in 

the lowest order with respect to ^i, ^2 and 9 

/ N ^9 ^ ^ exp(2z(r — 2r„)) 
< 0<r)>.^ 9^ ,^^);;;)3 + ^-^- r>r., (17) 

One can see that the amplitude of the echo is maximal at the time r = 2Tp. Note that 
the echo amplitude given by Eq.(p!7D is proportional to the product ^1^2 in contrast to the 
dependence < 7 >e~ typical for the phonon echo mechanism (for the small ^1, ^2) 
as a response to the external resonant forces |p. However, the spin echo type mechanism 
described above also leads to the terms ~ omitted from Eq. (p^) . 

According to the conventional model, the phonon echo originates from the change of 
the rotation frequency of the A-vectors due to the external pulses. This mechanism is 
also present in our model if we take into account the change of the energy E in each time 
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interval. In this case, the main contribution comes from the second term on the r.h.s. of the 
Eg. (p!5D ■ However, this effect is proportional to C,i^2 ^^"^ vanishes for small pulse intensities. 
Analytical calculations in this case become rather cumbersome and we do not present them 
here. Since we disregard the change of the rotating frequency due to the second pulse, the 
approximation leading to Eq.(PTD is valid if AE2Tp ^ C,2EoTp << 1. This condition can only 
be satisfied in the limit ^2 << 1 for large evolution time (rp >> 1). This is a consequence 



of the asymptotic nature of the proposed analytical solutions []T4 



Note, that the echo effect analyzed above is a stimulated classical revival of breathing 



oscillations in contrast to the effect of quantum revival caused by the discrete nature of 
the energy levels in quantum systems. We estimated a characteristic time of the quantum 



revival (or = ootr) by extending the method of Ref. ||T2[ to finite temperatures. The 
ratio Tc/Tr is AhuQ/SnT and therefore » t^. in the classical limit. Since the classical echo 
mechanism described above takes place in the time domain r 2rp, the quantum revival 
will not contribute to the effect if Tp < r,,. For the value of = 0.03 used in our numerical 
simulations and assuming that the number of levels in the potential well (|I|) is ~ 100, we 
obtain ~ 30 and ~ 300. Since 2Tp ^ 200, the above mentioned condition is satisfied. 



IV. CONCLUSION 

We investigated the dephasing and stimulated revival (echo) of the breathing oscillations 
in optical lattices subject to two-pulse parametric excitation created by the abrupt change 
of off-resonant dipole optical potential induced by the laser field. The dephasing originates 
from the energy dependent corrections to the oscillation frequency and due to the thermal 
distribution of energy. The shape of the echo response was obtained by numerical simulation 
using the equation of motion for the atoms with subsequent averaging over the initial thermal 
distribution. An approximate analytical analysis demonstrates that the suggested echo 
mechanism incorporates the essential physics features of spin(photon) echo where echo effect 
is caused by the change of the phase of the oscillations due to the second pulse, as well as the 
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phonon echo where the effect is caused by the change of the oscillation frequency. An optimal 
control analysis would be desirable (e.g. as proposed for the spin systems [|T5|) for optimizing 
the echo response in view of the rapid development of novel experimental techniques (e.g. 
[0)for the coherent manipulation of the atomic motion in optical lattices. 
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FIGURES 

FIG. 1. Fig.l. Numerical simulation of the echo effect discussed in the text. The parameters 
for the simulations are Tp = 105, (^i = 0.2, ^2 = 0.1, = 0.03. 

FIG. 2. Fig. 2. The maximum echo amplitude as a function of the intensity ^2 of the second 
pulse. 

FIG. 3. Fig.3. Illustration of the echo mechanism. Complex vectors Ai and A[ correspond to 
the different values of the atom initial energy. For < r < r^, the vectors rotate counterclockwise 
with the frcqTicncies ui'i > toi and accumulate the phase shift. At the moment t = Tp, the second 
pulse is applied to the system changing the oscillation amplitude. The new amplitude A2 is a 
supperposition of Ai and its complex conjugate Ai* For t > Tp, the phase difference between 
vectors Ai and Ai decreases leading to the echo effect. 
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